ITEP-31/98 
August 1998 



Kink-antikink interactions in the double 

66" 

sine-Gordon equation and the problem of 
a,: resonance frequencies 

; V. A. GANP 



Moscow State Engineering Physics Institute (Technical University) , 
Kashirskoe shosse, 31, Moscow, 115409, Russia 



CO 

> 
in 

O o,nd 
G\ • 

O 1 Institute of Theoretical and Experimental Physics, Russia 

00 ' 
On 

|; A. E. KUDRYAVTSEV f 

^ ; Institute of Theoretical and Experimental Physics, 

o 

o 



B.Cheremushkinskaja, 25, Moscow, 117259, Russia 



^ ■ Abstract 



We studied the kink-antikink collision process for the "double sine-Gordon" (DSG) equation 
in 1+1 dimensions at different values of the potential parameter R > 0. For small values of R we 
discuss the problem of resonance frequencies. We give qualitative explanation of the frequency 
shift in comparison with the frequency of the discrete level in the potential well of isolated kink. 
We show that in this region of the parameter R the effective long-range interaction between 
kink and antikink takes place. 
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1. Introduction 

Resonant energy exchange mechanism, that we shall consider in this our paper, was orig- 
inally observed in the kink-antikink collisions for the A02-theory. To examine such process 
one should consider initial configuration in the form of kink (K) and antikink (K) placed at 
x = ±xo (xo ^> 1) moving toward each other with some velocities Vi. It was found that there is 
critical value of the initial velocity v cr = 0.2598 and at i>j > v cr inelastic KK scattering takes 
place while at v-i < v cr kink and antikink form a bound state. This bound state then decays 
into small oscillations [1]. 

Later on, when the KK collision process was studied more careful, so-called escape windows 
were found [2] in the range of the initial velocities < v cr . Escape windows are nothing more 
than some values of the initial velocity Vi = v n at which kinks escape to infinity after second 
collision instead of forming a bound state. This phenomenon was semiquantitative explained in 
Ref. [2]. The point is that the A02 _ theory kink excitation spectrum has one zero (translational) 
and one non-zero (shape) mode with the frequency lo\ = \Js/2. It was observed that the 
following condition is satisfied with a reasonable accuracy: 

w 1 T 12 (v n ) = 6 + 2Tm, (1) 

where T 12 is the time interval between the two collisions of the kinks, n is integer, 5 is some 
constant phase. During the first KK collision a part of their kinetic energy is transferred to 
excitation of the kink discrete mode uj\. Therefore, kinks cannot escape to infinity and only 
go away at some distance and collide again. If condition (1) is satisfied, part of energy which 
conserved in the mode u 1 , is returned back to kinks translational mode (kinetic energy) and 
kinks can overcome the mutual attraction and go to infinity. Just this phenomenon was named 
" resonant energy transfer mechanism" . 

Note, that " higher orders" escape windows were also found. In these cases K and K escape 
to infinity after three or more collisions. For more detailed information about solitary wave 
interactions in the classical field theory see Review [3]. 

Let's now turn our attention to the system which we shall investigate in the present paper. 
The double sine-Gordon (DSG) equation can be obtained from the Lagrangian of the form 

with the potential 



^(0) = - 1+4 ^| ^cos0-cos-j . (3) 

Parameter 77 may be assigned any arbitrary real value (—00 < rj < +00). From Lagrangian (2) 
we get for the real scalar field (f>(x,t) in (1+1) dimensions the following equation: 

d 2 d) d 2 d) 2 / <A 
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In the present work we shall consider the range rj > 0. In this case it is suitable to introduce 
parameter R related with r] by the equality: 

n — — sinh 2 R . 
4 

Equation (4) has static topological solution in the form of 47r-kink (antikink): 

sinh x 

4>K(K)( X ) — 47rn ± 4 arctan — - . (5) 

v ; cosh R 

The sign "+" corresponds to the case of kink, "— " - to the case of antikink, n is integer. Eq. (5) 
can be rewritten in the form: 

4>K{R){x) = 4?™ ± [(f) S GK(x + R) - 4>sgk(R - x)\ , (5a) 

where 4>sgk{%) = 4 arctan exp(x) is the sine-Gordon (SG) equation 27r-soliton. From Eq. (5a) 
the physical meaning of the parameter R becomes clear: DSG kink can be interpreted as a 
superposition of two SG solitons, separated by the distance 2R. 

The KK collision process at a variety of values of the parameter 77 and the initial velocity 
Vi was studied in details in Ref. [4]. As for the A02-theory case, there is some critical velocity 
v cr below which kinks form a bound state decaying into small oscillations. Note, that in the 
DSG case this critical velocity is a function of the parameter (rj or R) [4]. 

It was found that in the DSG system the resonant energy exchange mechanism also takes 
place. As a consequence there is a system of escape windows at some values of rj. Note, that 
there is one important difference in kinks collision processes between A02 an d DSG models. In 
the first case kinks cannot pass through each other at a valuable distance, while in the second 
they can travel to infinity after passing through each other. This difference is a consequence of 
the different structure of the potential V(<f>). 

In Ref. [4] different values of the parameter R were studied. At it! = 1.2 a typical picture of 
escape windows was similar to the A02-theory case. However, at smaller R, namely at R = 0.5, 
a new phenomenon was observed in the KK collisions - so-called quasiresonances. The essence 
of the phenomenon is in the following. At all velocities Vi < v cr we get capture and formation 
of KK bound state. But the time between the second and third collisions T 23 as a function of 
the initial velocity Vi has a series of well-defined maxima, see Fig. 1. Such behavior of T 2 2,(vi) 
means that the resonant energy exchange mechanism appears in the system, but at the same 
time the energy returned to the translational mode during the second collision is not enough for 
kinks to escape to infinity after the second collision. Besides, it turned out that the frequency 
of oscillations in which a part of kinetic energy is transferred while kinks passing through each 
other for the first time, i.e. oo\ in Eq. (1), is smaller than the frequency of the localized DSG 
kink excitations at given R. 
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2. General approach 



Up to now, investigating the resonant energy exchange mechanism, people involved the 
localized excitations over an isolated kink (or antikink). As it will be shown, in some cases 
of the KK scattering such approximation is not valid. In these cases one should consider the 
spectrum of small excitations for the KK system as a whole. For this purpose let's look for 
the solution of Eq. (4) for the field in the form: 

(f)(x, t) = 2ir + 4>k{x - x ) + 4>k{ x + x o) + 5(f)(x, t) . 

Such configuration corresponds to the kink and antikink placed at x — ±x plus some small 
perturbation 8(f)(x,t), \5<f>\ <C 1. Taking into account that 0x and are solutions of Eq. (4), 
we get for S(f> the following linearized equation: 



+ 50 



d 2 V 



d<p 2 



= Q(x,x ) , 



where 



dV 

Q(x,x ) = — 



+ 



dV_ 

90 



dV_ 



(6) 



(7) 



<f>=2TT + (f>K+(f>K 

The explicit form of the function Q(x,xo) is rather cumbersome (see Appendix), but never- 
theless we can make several general notes. Inhomogeneity Q(x,x ) in Eq. (6) is a consequence 
of the fact that the configuration "kink+ antikink" is not a solution of Eq. (4). The function 
Q(x,xo ) characterizes overlapping of the kink and antikink, because <Pk(x — xq) and (p^(x + x ) 
are exact solutions of Eq. (4) when taken separately. Obviously, Q(x,x ) is an even function 
of x and x and it falls down exponentially when xo increases. At fixed xo as a function of x 
Q(x,x ) looks like two bumps with maxima at x — ±x . 

Let's now find the excitation spectrum for 5 (p. For this purpose we take Eq. (6) with zero 
right-hand side and look for S(f> in the form: 

S(j>{x,t) = e^xix) ■ 

Then for the function x{ x ) we get the following differential equation of the Schrodinger type: 

-x" + U(x,x )x = uj 2 x , (8) 



where 



U(x,x ) = 



d 2 V 



d<p 2 



(9) 



<t>=2n+4> K +(j> R 

The explicit form of the potential U{x,xq) is rather complicated (see Appendix) and depends 
crucially on xq. Note, that the shape of this potential depends on the parameter (R or 77) 
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and U(x,x ) — > 1 when x — > ±00. Hence, u> < 1 form the discrete excitation spectrum, and 
c<j > 1 - the continuum one. In the limit x ^> 1 C/ (x, x ) as a function of x looks like 
two identical potential wells, separated by the distance 2x . Each well contains one or more 
discrete levels which correspond to the localized excitations of the solitary kink (antikink). In 
the collision process DSG kinks pass through each other, i.e. Xq decreases to zero and then 
starts to increase again. At small xq the distance between the wells is small and the discrete 
levels are not independent. With kinks moving toward each other from the infinity the levels 
begin to split and then at x ^ 1 mutual potential of the system KK is quite different from 
one of the solitary kink (antikink). 

It is worth mentioning that taking into account of both wells is also necessary in cases 
when in each potential well there is a discrete level with small binding energy situated near the 
continuum. In such cases one should take into account overlapping of the wave functions in 
both wells even at xq ^> 1. It means that under some conditions long-range interaction between 
kink and antikink appears in the system. 

In what follows we will show, that within such approach the phenomenon of quasiresonances 
observed in the DSG system at R = 0.5 in Ref. [4] may be simply explained. We will also argue 
that the cause of the quasiresonances is just the resonant energy exchange mechanism, that 
leads to escape windows at some other values of R. Moreover, there is some intermediate region 
of R where quasiresonances and escape windows appear together. 

3. Small R 

In Ref. [4] quasiresonances were observed at it! = 0.5. We performed similar calculations 
and obtained analogous curve T^Uj), see Fig. 1. Besides that, we have investigated the KK 
collision process at R — 0.4 and R = 0.6. At R — 0.4 (Fig. 2) we get a picture of quasiresonance 
peaks analogous to the case R = 0.5. At R — 0.6 (Fig. 3) there seem to exist escape windows in 
places of some peaks on the curve T2s(vi). It confirms that quasiresonances and escape windows 
are phenomena of the same nature and with the parameter R increasing some quasiresonance 
peaks transform into the escape windows. At some intermediate values of R both phenomena 
are presented, and in further increasing of R only escape windows survive. At R — 1.2 in 
Ref. [4] a perfect picture of escape windows and no quasiresonances were observed. 

To answer the question why at given R quasiresonances or escape windows appear it is 
required, generally speaking, to solve Eq. (6) for 5(f) with the right-hand side. At the same 
time we can suggest some true-like hypothesis. Each bump of the source Q(x,x ) is localized 
on size of order of 1 (see Appendix). In the case of small R ~ 0.4 — 0.6 the first excited level 
in the well is not well-localized (binding energy is small). Therefore, integral of overlapping of 
Q(x, xo) and the wave function of the excited state is small. It corresponds to the fact that the 
part of the kinetic energy transferred to the discrete mode lo\ is small, and hence loss of energy 
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due to radiation is large. The situation changes with increasing of R. The binding energy of 
the first excited level is increasing, and for R — 1.2 the first excited level in the well is already 
well-localized. Because of this reason the character size of the wave function is of the same 
order as the source one. Hence the energy transfer mechanism is more effective in this case. 

From analysis of the quasiresonance peaks of the T 23 (fj) plot for R = 0.5 (Fig. 1) it follows 
that the frequency of the discrete mode being excited is approximately equal to £o\ = 0.945. 
At the same time in the well corresponding to one kink excitations there is discrete level with 
frequency u\ = 0.967. As it will be shown, this deviation is not incidental and may be easily 
interpreted within our approach. 

At small R the system is close to the pure sine-Gordon. Therefore, the critical velocity v cr 
is small (v cr = corresponds to the pure sine-Gordon) and potential (9) in the Schrodinger 
equation (8) has one discrete level situated near the continuum (in the pure sine-Gordon case 
there is only zero mode). Presence of a shallow level implies that the corresponding wave 
function falls down slowly with the distance from the well. Due to this reason while studying 
the KK collision process it is necessary to take into account the fact that the wells affect each 
other even at large distance. It leads to changes in the excitation spectrum. In Fig. 4 we show 
how the excitation frequency of the KK system depends on the distance between K and K 
(this distance is equal to 2x ). From the plot it is seen that even at x 3> 1 there present some 
visible splitting of the higher discrete level. In presence of the second well this level oj\ = 0.967 
splits into two sublevels: the higher with £j° dd > 0.967 and the lower with Cuf ven < 0.967. In 
the collision process the lower one is excited because the corresponding wave function is even. 
Moreover, the higher level may disappear during the collision. From Fig. 4 one can see, that 
this higher level appears from the continuum uj > 1 at some critical distance 2(x )cr 1 
between kink and antikink. So, if the original level in a single well lies near the continuum (as 
it happens at small R), then (x ) cr is very large. At the same time the level with ojf ven exists in 
a relatively wide interval of distances xo between kink and antikink. Namely this level is excited 
during the KK collision because of the resonant energy exchange mechanism. In this case the 
frequency cu 1 in expression (1) is indeed the averaged over different x the frequency Cj{ ven . It 
is smaller than the frequency of the discrete mode of isolated kink, what is in correspondence 
with the numerical simulations. 

4. Numerical calculations 

We solved the second order partial differential equation (4) numerically on the lattice with 
Ax = 0.01. Initial conditions were taken in the form of kink and antikink (5) situated at 
x = ±20 moving towards each other with velocities =Ft> j respectively. Moments of kink and 
antikink passing through each other were fixed via field behavior at the origin x — 0. 

To find discrete levels in the potential (9) we used the fact that the wave function falls 
down exponentially at large distances. We took solution of the Schrodinger equation in the 
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form x ~ exp (xy/l — uj 2 ) at x — —50 and solved numerically stationary equation (8). As a 
result we got x at x = 50 as a function of uj. When uj does not correspond to the discrete 
level, x grows exponentially with x at positive x's. But if uj coincides with a discrete level of 
the potential, then x is exponentially suppressed at large x. In real computations we observed 
that x( x — 50) changed its sign when uj passed a discrete level. Note, that this method being 
applied to searching of a shallow level does not return a good result. In such a case one should 
take more distant starting and ending points. The origin of the problem is in the following: for 
a shallow level uj — > 1, and exponents fall and grow very slowly with increasing of x. 

Conclusion 

This our paper presents qualitative and semiquantitative explanation of the phenomenon of 
quasiresonances in collisions of kink and antikink of the double sine-Gordon equation at small 
R. It is shown that the resonant energy exchange mechanism being applied in its previous form 
gives not satisfactory results for frequencies. 

It was shown that the resonant energy exchange between kinks' translational mode and the 
discrete excitations of the KK system as a whole takes place. At small R it is essential because 
of long-range interaction in the system caused by the presence of a shallow level in the discrete 
spectrum of excitations of an isolated kink (antikink). 

The proposed mechanism explains qualitatively the decrease of the resonance frequency uj\ 
in Eq. (1) at small R in comparison with the discrete frequency of an isolated kink. 
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Appendix 



In equation 



d 2 



d 2 4> dV 

+ — = 0, 



dt 2 dx 2 ' d<j) " ' 

where V(0) has the form (3) let us substitute 4> = 2n + 4> K + 4>r + S4>- Taking into account 
that 4>k and (fig are exact solutions of (Al) and linearizing with respect to 6(f) we get: 



d 2 5(j) d 2 



d 2 V 



dt 2 dx 2 



5(fi = Q(x,x ) , 



(A2) 



2n+4>K+4>K 



where Q(x, Xo) is given by Eq. (7). If we substitute explicit expressions for kink and antikink 
situated at x — ±x respectively, then we get: 



Q(x,x ) 



+ 4:7] 



- 4) 
(1 + s 2 ) 2 



1 + 4.7] L(i + 4)(i + s 2 

2xA s + (l-s 2 + ] 



.1 + 4 



here 



s± 



(i + 4) 2 

sinh (x ± Xo) 



1-4 

1 + s 2 



(A3) 



coshi? 

[j] and R are related by r] = (1/4) sinh 2 i?.) 

Let us use Q(x,xq) = in (A2) and substitute 5(f)(x,t) = x( x ) ex P (^0> then we obtain: 



where 



U(x,x ) 



d 2 V 



Insert here explicit expressions for kink and antikink. As a result we have 



U(x,x ) 



1 - 4 1 - s 2 



+ 



+8r/ 



i + s 2 _i + 4 (i + s 2 )(i + 4) 

2 



1 +4r] 

l - 4 i - 4 
i + s 2 _i + 4 + (i + s 2 _)(i + 4) 



4s s. 



- 4r] 



(A4) 
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Figure captions 



Fig. 1. The time T 12 between the first two KK collisions (dashed curve) and the time T 23 
between their second and third collisions (solid curve) as functions of the initial velocity V{ for 
R = 0.5. 

Fig. 2. The time T 12 between the first two KK collisions (dashed curve) and the time T 23 
between their second and third collisions (solid curve) as functions of the initial velocity Vi for 
R = 0.4. 

Fig. 3. The time T 12 between the first two KK collisions (dashed curve) and the time T 23 
between their second and third collisions (solid curve) as functions of the initial velocity v i for 
R = 0.6. Arrows denote probable positions of the escape windows. 

Fig. 4. The excitation frequency as a function of the initial half-distance between K and 
K. R = 0.5. 
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